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ABSTRACT 

This is the third in a series of papers that present observations and resuhs for a sample 
of 76 ultra-steep-spectrum radio sources designed to find galaxies at high redshift. Here we 
present multi-frequency radio observations, from the Australia Telescope Compact Array, for 
a subset of 37 galaxies from the sample. Matched resolution observations at 2.3, 4.8 and 
6.2 GHz are presented for all galaxies, with the z < 2 galaxies additionally observed at 
8.6 and 18 GHz. New angular size constraints are reported for 19 sources based on high 
resolution 4.8 and 6.2 GHz observations. Functional forms for the rest-frame spectral energy 
distributions are derived: 89% of the sample is well characterised by a single power law, whilst 
the remaining 11 % show some flattening toward higher frequencies: not one source shows any 
evidence for high frequency steepening. We discuss the implications of this result in light of 
the empirical correlation between redshift and spectral index seen in flux limited samples of 
radio galaxies. Finally, a new physical mechanism to explain the redshift - spectral index 
correlation is posited: extremely steep spectrum radio galaxies in the local universe usually 
reside at the centres of rich galaxy clusters. We argue that if a higher fraction of radio galaxies, 
as a function of redshift, are located in environments with densities similar to nearby rich 
clusters, then this could be a natural interpretation for the correlation. We briefly outline our 
plans to pursue this line of investigation. 

Key words: surveys - radio continuum: general - radio continuum: galaxies - galaxies: 
active 



1 INTRODUCTION index selection. 



Until the late 1970s the search for high redshift radio galaxies 
(HzRGs) was restricted to optical spectroscopy at the position of a 
radio source whose optical host galaxy is undetected on, or at the 
detection limit of, the Palomar Observatory Sky Survey (POSS) 
plates (e.g .Kristianetal. 19 74. 1978; Smith & Spinrad mOi 
[Gunn et al]ll98ll: Laing et alJl l983: Soinrad & Diorgovs]dl ll984l: 
[Perrvman e^i'^^84nSpin]radet aljl985tV Indeed, the fir st galaxy 
disco vered above redshift one was found in this way iSpinrad 
Il982h . 

The field was revolutionised when it was realised that the 
radio sources without optical counterparts on the POSS plates 
tended to have the steepest spectral indices q (where oc v") 
measured between 178 a nd 1415 MHz iTielens et alj 1 19791 : 
iBlumenthal & MilevI [19791) . Since then most radio selected sam- 
ples designed to find HzRGs have used some degree of spectral 



Three explanations have been put forward in the literature to 
justify the trend for HzRGs to have steeper spectral indices; the 
so-called z — a correlation. The first is based on the observation 
that the spectral energy distributions (SEDs) of radio galaxies 
typically steepen toward higher frequencies. For a fixed set of 
observing frequencies, a radio galaxy at higher redshift will be 
sampled at higher rest-frame frequencies where, on average, it will 
exhibit progressively steeper spectral indices: a k-correction. It has 
also been suggested that the SEDs of HzRGs will be intrinsically 
steeper thanks to enhanced inverse Compton (IC) losses of the 
relativistic electron population against the cosmic microwave 
background (CMB) p hotons w hose energy density increases as 
il + z)^i Ki-olik & C herJl99lh . These two explanations have been 
adopted in the literature over recent years as the main contributors 
to the correlation ( De Breuck et al. 20.00. 2002; Pedani 20031 
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ICohen et a l. 2004: J arvis et alJl 2004Y The third explanation holds 
that the z — a correlation itself is indirect, reflecting an intrinsic 
correlation between ra dio luminosity and sp ectral index coupled 
to a Malmquist bias iSlundell et alj Il99 9h. Regardless of the 
underlying reason for the correlation, the steep-spectrum technique 
has enjoyed tremendous success over the years, with at least 30 
radio galaxies above 2 = 3 found to date using some form of 
spectral index culling. 

I n the first paper of this series (paper I; IPe Breuck et alJ 
l2004h . we defined a sample of 76 high redshift radio galaxy 
candidates selected between 843 and 1400 MHz on the basis of 
their ultra steep radio spectral index (USS; a < —1.3). In the 
second paper of the series (paper II; IPe Breuck et aljEzOOd) we 
reported the results of optical spectroscopy for 52 sources in the 
sample which included 35 spectroscopic redshifts. The aim of this 
paper III is to better understand the radio spectral characteristics of 
our sample of USS radio galaxies, over a broad range of redshifts. 
To this end, we report multi-frequency radio observations from 
the Australia Telescope Compact Array (ATCA) and present 
rest-frame radio SEDs for 37 galaxies in the sample. 

This paper is organised as follows. Some relevant synchrotron 
physics is given in O The subset of SUMSS-NVSS USS ra- 
dio galaxies we observed for this paper is defined in Js] Multi- 
frequency ATCA observations and data analysis follow in ij4| The 
results of the observational program are presented in [js] which in- 
cludes rest-frame radio SEDs. The remainder of the paper is dis- 
cussion based. S|6|and focus on the implications of our results 
for existing interpretations of the redshift - spectral index correla- 
tion. Penultimately, in S|8| we speculate on an alternative physical 
mechanism for this correlation and concluding remarks are given in 
S|9| Throughout this paper, we assume a flat A-domi nated cosmol- 
ogy w ith Ho = 71 kms^^Mpc"^ and = 0.73 iSoergel et alj 



in vacuo: 

P P^du 

(3) 

" Gttc ^" ' c ' 
where a is the acceleration and /io is the permeability of free space. 
For circular motion where a is the angle between the velocity and 
acceleration vectors: 

carB'^'y^ . 2 

P = ^sm a (4) 

MO 

where ar is the Thomson scattering cross-section for electrons. 
The observed frequency spectrum of synchrotron emission is often 
characterized by a power law. This spectrum is interpreted as origi- 
nating from a distribution of relativistic electrons with a power law 
energy distribution of the form: 

N{E)dE = NoE^dE (5) 

where F is a constant power law index and A'^o is the total elec- 
tron number density. Substituting Equations|4|and|5|into|2| and as- 
suming an isotropic distribution of electrons so that (sin^a) = |, 
the power per unit frequency per unit solid angle averaged over an 
emitting volume V is given by: 

= -— car — — NqE' V , (6) 

DTT ZflQ V 

with the relationship between 7, B and v being given by Equa- 
tionQ Also, since E — 'ynioC?, then from EquationQwe can see 
that E^ (X . Therefore, the radiated power at a specified observ- 
ing frequency obeys the proportionality: 

P - r + 1 

P„ (xu , a = — - — (7) 

where we define a to be the spectral index of the synchrotron radi- 
ation. 



2 SYNCHROTRON RADIATION 

A relativistic electron spiralling around a magnetic field line emits 
synchrotron radiation expressed as: 



eB 
2Tvmo 



(1) 



Equation Q which relates the specific frequency u of the syn- 
chrotron radiation for relativistic particles with rest mass mo in a 
magnetic field B with a Lorentz factor 7, is the fundamental re- 
lationship for synchrotron radiation processes in astrophysics. The 
emission from an ensemble of electrons with energies between E 
and E + dE can be approximated by a delta function at a frequency 
given by Equation^ The synchrotron emissivity is: 



= J S{x~ E')P^N{E) dEdx 
= P„N{E')dE' 



(2) 



where P„ is the power per unit frequency emitted by one electron 
from an ensemble of A'^ electrons per unit volume per unit solid 
angle with energies in the range from E' to E' + dE' . Pi, can be 
determined from the relativistic generalization of the Larmor for- 
mula which describes the total power P radiated by a point charge 



2.1 Electron Energy Losses 

We now consider energy losses to the initial electron energy dis- 
tribution, which lead to a departure from Equation |5| The general 
form of the energy loss f{E) can be represented by: 



cp{E) 



-C-riE^ ^E' 



(8) 



where ( represents the electron energy losses which are essentially 
independent of E, rj represents the losses proportional to E and ^ 
the losses proportional to E^ . We consider synchrotron and inverse 
Compton cooling which contribute to ^ and adiabatic losses which 
contribute to rj. 

2.1.1 Synchrotron Losses: the ^ term 

In a constant ambient magnetic field the radiating electrons will 
lose energy at a rate which is given by Equation |4| The timescale 



for an electron with initial energy E2 
its final state Ei = 71 mo (? is then 



72 mo c to lose energy to 



1 



dE 



c (tt B"^ 71 



72 ' 



(9) 
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In a magnetic field, this electron will radiate at continuously lower 
frequencies over this timescale (Equation 0. For a continuous in- 
jection model I Kardashev 1962; Carilli et al. 1991 ), the initial elec- 
tron energy distribution will steepen from — F to — (F + 1) above 
a critical, or 'break' energy. This will result in the associated syn- 
chrotron SED remaining a power law but with a spectral steepening 
from a to a — i above a critical frequency which continuously de- 
creases with increasing time. In the limiting case where there was a 
single injection of electrons, then the SED steepens exponentially 
when t > T. 



2.1.2 Inverse Compton Scattering: the ^ term 

If the energy density of the radiation field is high enough, the syn- 
chrotron emitting electrons will collide with the ambient photons, 
suffering energy losses as they impart their momentum to the radia- 
tion field. As for synchrotron losses, the rate of energy loss is again 
proportional to the square of the electron energy, meaning that the 
resultant synchrotron spectral index will again steepen from a to 
Q — ^ . The ambient photon field consists of the synchrotron pho- 
tons themselves (in which case the process is called synchrotron 
self-Compton, SSC), the cosmic microwave background (CMB) 
photons, and the infra-red (IR) photons from the AGN and galaxy. 
Thus, considering inverse Compton cooling as well, Equation|9|be- 
comes: 



El 



dE 



E2 



COT (S2+B„(^)2) J^^ £2 

MO 

CGT (B2+B„(z)2) ^-yi 



(10) 



where Bic (z) is the equivalent inverse Compton magnetic field at 
redshift z. The ratio of the magnetic energy density 

Ub - - — 
to the radiation energy density: 



(11) 



by some linear expansion factor I. In the absence of tangled mag- 
netic fields, conservation of magnetic flux (the third adiabatic in- 
variant) gives. 



(14) 



4.1 X 10-^\l + z)^ + Ussc + (7,R, 



(12) 



From Eauationsm il3l and ll4l the synchrotron SED will be shifted 
to lower frequencies by a factor of In addition, the radiated 
power of an electron, as given by Equation |4| is proportional to 
7^ B^. If the frequency of this radiation was fixed, then the am- 
plitude of the radiated power would decrease by l^, because B^ 
decreases by and 7^ decreases by l^. However, since the fre- 
quency of the radiation decreases as (from EquationQwhere B 
decreases by l'^ and 7^ still decreases by l^), then the actual radi- 
ated power at a fixed emission frequency decreases by l^. This is 
a somewhat over-simplified model which we use to make the point 
that the flux density, at a given frequency, from the lobes of a radio 
galaxy will decrease over time at a rate which is governed by the 
pressure of the environment they are expanding into. Thus, a ra- 
dio galaxy selected from a given flux-limited survey, whose lobes 
remain pressure-confined by a relatively high ambient gas density, 
may fall below the survey flux limit if it was, instead, subjected 
to larger adiabatic losses expanding into a more rarefied environ- 
ment. For more comprehensive and realistic models of classical ra- 
dio galaxies, we refer the reader to Blundell et al. 1 1999.) and refer- 
ences therein. 



3 TARGET SELECTION 

We selected 37 (of the 76) SUMSS-NVSS USS radio galaxies for 
follow-up multi-frequency radio observations. This subset consists 
of the 28 galaxies whose redshifts had been spectroscopically con- 
firmed by 2004 September, and another nine galaxies without con- 
firmed redshifts. Seven of the latter group have extremely faint host 
galaxies (K > 19.88) and are likely to be at red shifts z > i ac- 
cording to the K — z rel ation for radio galaxies jPe Breuck et alJ 
I2OO2; Willott et 4||2003). The final two sources without redshifts 
were included in the subset purely due to their favourable hour an- 
gle during the ATCA observations when extra observing time b e- 
came available. The K ~ z relation from lDe Breuck et all2002l is 
used to estimate the redshifts of the nine galaxies without spectro- 
scopically confirmed redshifts. 



in units of erg cm~'^, will determine whether synchrotron losses 
or inverse Compton losses dominate the electron energy loss spec- 
trum. 

2.1.3 Adiabatic Expansion Losses: the rj term 

We follow closely the derivation o f IScheuer & WilliamsliT963) . An 
initial energy distribution of electrons described by N{Ei) in a vol- 
ume of radius Ri, undergoing adiabatic expansion to a volume of 
radius Rf is governed by the invariance of TV^~^, where T is 
the temperature (energy), V is the volume and F = | is the usual 
adiabatic index for relativistic gas. Then, 

N{Ef)=N{E,)i^f^^-'\ (13) 

where N{Ef) is the energy distribution after expansion. Therefore, 
under adiabatic expansion, N{Ef) will decrease relative to N{Ei) 



3. 1 Revised 843 MHz flux densities 

Whilst determining the SEDs which are presented in il5.2l we 
noticed that almost all flux densities were well defined by a 
single pow er law except for th ose at 843 MHz from the SUMSS 
catalogue ilMauch etalj|2003h which were systematically larger 
than expected. Re-measurement of these fluxes using the SUMSS 
images showed that the catalogued values were high by a few 
percent. 

The reason for this over-estimation of the flux density is now 
understood. It only occurs at the very northern edge of the SUMSS 
catalogue, which is precisely where our USS sample has been 
selected, and arises from an incorrect assumption inherent in the 
source fitting algorithm (VSAD) which is optimised for a circular 
beam. This error is exacerbated for our sample because the USS 
selection favours sources with the largest flux over-estimate at 
843MHz. On average, the over-estimation is at the 4% level and 
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will be corrected in the next release of the SUMSS catalogue. 

We have manually inspected the relevant SUMSS mosaics and 
revised the 843 MHz flux densities for the 37 sources considered 
in this paper. The revised flux densities are listed in Column 5 of 
Table 2. A consequence of these adjustments is that several of the 
original USS sources now have spectral indices flatter than the 
a843° = -1-3 USS threshold. 



4 OBSERVATIONS AND DATA REDUCTION 

We used the ATCA to measure accurate radio positions and mor- 
phologies for the 16 highest redshift candidates in the sample: these 
included galaxies with (i) spectroscopic redshifts z > 2, (ii) red- 
shifts predicted from the K ~ z relation. These sources are marked 
with an asterisk (*) in Column 1 of Table 2. In order to measure 
total flux densities for all 37 sources in our sample, we further 
obtained matched angular resolution observations in the 2.4 GHz, 

5 GHz, 8.6 GHz and 18 GHz bands; a journal of the radio ob- 
servations is given in Table Q In all cases, we utilised the stan- 
dard continuum correlator configuration which has dual frequency 
mode with two independent 128 MHz bands (frequency I and 2 
in Table 0. For all frequencies, the flux density scale was deter- 
mined by observations of the standard ATCA primary calibrator 
PKS BI934-638. We observed PKS BI921-293 to correct for the 
instrumental bandpass, and a suite^ of compact bright active galac- 
tic nuclei (AGN) located in close proximity to the science targets to 
track variations in the atmospheric phase during the observations. 



4.1 A note on observing strategy 

We observed our sources using multiple short observations, making 
between four and eight snapshots per source spread over a range of 
hour angles. With such a large number of targets, telescope slew- 
ing time becomes appreciable so that some thought was given to 
the trade-off between optimising the uv-coverage and maximising 
the on-source time. The relatively large field of view at 2.3 GHz, 
and to a lesser extent at 4.8 GHz, results in many confusing sources 
within the primary beam. At these frequencies, the flux densities of 
the targets are still relatively bright and our experience has shown 
that it is optimal to maximise the uv-coverage to minimise effects 
of confusion. Therefore, for frequencies of 4.8 GHz and below, tar- 
gets were observed for three minutes each and we observed a phase 
calibrator every 45 minutes. On the other hand, the number of con- 
fusing sources in the primary beam becomes negligible at frequen- 
cies above 6.2 GHz, but the target sources are now extremely faint 
and the atmospheric phases less stable. For these observations, we 
sacrificed uv-coverage for sensitivity by making fewer snapshots 
with longer integration times (10-15 minutes) and we observed a 
phase calibrator every 30 minutes. 



frequency (RF) interference were removed; this was mostly a 
problem at 2.496 GHz, and to a lesser extent at 2.368 GHz. At 
17.856 GHz and 18.496 GHz, we removed data when the rms 
path length through the atmosphere (measured from a two-element 
interferometer tracking a geo-stationary satellite) was sufficient 
to cause decorrelation during the individual scans (approximately 
15% of the time). The bandpass was determined along with the 
absolute flux density scale, and both solutions were applied to the 
secondary calibrators whose complex gains and phases were deter- 
mined and in turn applied to the target sources. C ontinuum imaging 
was p erformed using multi-frequency synthesis jSault & Wieringal 
1 1 994ft prior to deconvolution. 



As mentioned in il4. II we were mostly limited by sidelobe 
confusion at frequencies below 6.2 GHz and limited by sensi- 
tivity above 6.2 GH z. We used images of the associated NVSS 
jCondon et ajl99 8h field obtained from its online image server to 
constrain the deconvolution algorithm to the 'real' sources in the 
field (as opposed to spurious sources due to phase errors or bright 
sidelobes), and accordingly restricted the deconvolution algorithm 
to tight regions surrounding each source. The sensitivity of the 
NVSS was sufficient to apply this constraint without biasing our 
results. This process is extremely time-consuming to do by-hand, 
and can also be difficult to reproduce systematically. As a result, an 
automated procedure has now been defined to select deconvolution 
(CLEAN) regions in ATCA images using the known positions of 
radio sources in the field. This is analogous to the FACES task 
which was independently implemented into the AIPS data analysis 
package several years ago. 

When the science target was unresolved, based on inspection 
of the deconvolved images and the visibilities, and when there were 
no confusing sources in the primary beam (usually at frequencies 
larger than 6.2 GHz), we determined the total flux density of the 
source from the peak of the 'dirty' map, which is more accurate 
than the deconvolved value. Otherwise, we measured total inte- 
grated flux densities from the deconvolved images. 



5 RESULTS 

Table 2 gives the results of these observations: 

Column I: NVSS source name. 
Column 2: Redshift from paper II. 

Column 3: 1.4 GHz radio luminosity L1.4, calculated using 



Li.. 



Wm- 



cm^ /ijy 



^10-^"^(l + 2) 



(15) 



where Dl (z) is the luminosity distance and 5i.4 is the flux density 
measured at an observing frequency of 1.4 GHz. 



4.2 Data reduction 

We emp loyed the M IRIAD radio interferometry data reduction 
package iSault et alJil995.) and used the same standard procedure 
to reduce all the observations. Initially, all data corrupted by radio 

1 PKS B 1954-388, PKS B2058^25, PKS B2211-388, PKS B2255-282, 
PKS BOOlO-401, PKS B0104-408, PKS B0153^10, PKS B0220-349 



Column 4: Spectral index between 843 MHz and 1.4 GHz. 
Column 5: Total flux density measured at 843 MHz. 
Column 6: Total flux density at 1.4 GHz from NVSS 
("Con don et alll998ft . 

Column 7: Total flux density measured at 2.4 GHz. 
Column 8: Total flux density measured at 4.8 GHz. 
Column 9: Total flux density measured at 6.2 GHz. 
Column IO:Total flux density measured at 8.6 GHz. 
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Table 1. Journal of the multi-frequency radio observations obtained with the ATCA 



Array 


Epoch 


Frequency 1 


Frequency 2 


Field of Viewt 


Beam Size^ 






GHz 


GHz 


arcmin 


arcsecxarcsec 


6A 


2004 December 10-12 


4.800 


6.208 


9.9 


3x6 


1.5A 


2005 April 16 


2.368 


2.496 


20.1 


17 X 32 


750A 


2005 April 29-30 


4.800 


6.208 


9.9 


17 X 31 


EW367 


2005 June 12-14 


8.568 


8.626 


5.5 


21 X 39 


H168 


2005 May 10-12 


17.856 


18.496 


2.7 


15 X 15 


fField of View (half power width of the primary beam) and Beam Size are quoted for Frequency 1. 



Column ll:Total flux density measured at 18 GHz. 
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Table 2. List of 37 sources in tliis sample in order of increasing Right Ascension. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Source 


z 


-^1.4ghz 


"843 


S843MHZ 


Si.4ghz 


S2.4GHZ 


S4.8CHZ 


S6.2GHZ 


S8.6GHZ 


SlSGHZ 






WHz~i 




mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


NVSS J001339-322445* 


0.2598±0.0003 


3.3x10^5 


-0.93 


248.4 ± 12 


155.2 ± 4.7 


82.5 ± 4.1 


47.2 ± 2.4 


37.4 ± 1.9 


26.1 ± 1.3 


10.7 ± 2.1 


NVSS J002 13 1-342225 


0.249±0.001 


3.8 xlO^* 


-1.43 


37.2 ± 2 


18.0 ± 1.0 


7.6 ± 1.2 


5.6 ± 0.6 


4.8 ± 0.5 


3.4 ± 0.4 




NVSS J002219-360728 


0.364±0.001 


7.8x10^* 


-1.34 


30.6 ± 2 


15.5 ± 0.7 


8.5 ± 0.4 


5.5 ± 0.3 


4.2 ± 0.3 


3.0 ± 0.3 


1.6 ± 0.3 


NVSS J002402-325253* 


2.043±0.002 


2.9x10^^ 


-1.34 


81.7 ± 4 


41.3 ± 1.3 


26.0 ± 2.0 


12.5 ± 0.6 


10.5 ± 0.6 






NVSS 1002627-323653* 


0.43±0.01 


3.8 xlO^-^ 


-1.25 


83.5 ± 4 


44.4 ± 1.7 


19.8 ± 1.0 


9.6 ± 0.8 


7.8 ± 0.5 


4.4 ± 0.4 


2.0 ± 0.4 


NVSS J011606-331241 


0.352±0.001 


1.1x10^^ 


-1.33 


43.3 ± 2 


22.1 ± 0.8 


9.0 ± 2.0 


3.6 ± 0.3 


2.3 ± 0.3 


1.6 ± 0.3 


0.9 ± 0.5 


NVSS JO 12904-3 248 15 


0.1802±0.0003 


3.8x102"* 


-1.39 


82 ± 4 


40.5 ± 1.7 


16.0 ± 4.0 


8.5 ± 2.0 


4.4 ± 1.0 


3.6 ± 0.5 


1.3 ± 0.3 


NVSS 1015232-333952* 


0.618±0.001 


4.6x10^'^ 


-1.05 


420.5 ± 21 


247.3 ± 8.7 


157.0 ± 8.0 


95 ± 5 


77 ± 4 


61 ± 3 


23.4 ± 4.7 


NVSS J015324-334117* 


0.1525±0.0004 


1.3x102* 


-1.21 


36.6 ± 2 


19.8 ± 0.8 


10.0 ± 2.0 


8.4 ± 1.0 


7.0 ± 1.0 


5.7 ± 0.6 


2.4 ± 0.5 


NVSS J015544-330633* 


1.048±0.002 


2.6x10^6 


-0.91 


61.8 ± 2 


39.0 ± 1.6 




13.0 ± 1.0 


10.0 ± 0.5 


7.7 ± 0.3 


3.2 ± 0.5 


NVSS J02 1308-322338 


3.976±0.001 


9.3 X 10^^ 


-1.48 


63.7 ± 3 


30.0 ± 1.0 


19.5 ± 1.5 


10.3 ± 1.0 


8.5 ± 0.6 






NVSS J030639-330432 


1.201±0.001 


3.9x10^6 


-1.69 


63.7 ± 3 


27.0 ± 0.9 


11.5 ± 1.0 


3.7 ± 0.3 


2.4 ± 0.3 


1.3 ± 0.2 


0.9 ± 0.5 


NVSS J202026-372823* 


1.431±0.001 


6.4x10^^ 


-1.27 


70.0 ± 4 


36.8 ± 1.2 




9.8 ± 0.5 


6.4 ± 0.3 


4.5 ± 0.3 


1.7 ± 0.6 


NVSS 1202140-373942* 


0.185±0.001 


2.0x102* 


-1.18 


35.4 ± 2 


19.5 ± 1.1 


11.0 ± 1.0 


6.0 ± 1.0 


5.6 ± 0.8 


4.6 ± 0.8 


2.2 ± 0.4 


NVSS 1202945-344812* 


1.497±0.002 


9.5x102'' 


-1.22 


97.5 ± 5 


52.5 ± 2.0 


32.0 ± 1.6 


17.1 ± 1.0 


13.4 ± 0.7 


11.2 ± 0.6 


3.5 ± 0.7 


NVSS J204420-3 34948* 


> 4t 


>7.7xl02'^* 


-1.52 


44.9 ± 2 


20.8 ± 0.8 


10.0 ± 0.5 


3.3 ± 0.3 


2.0 ± 0.3 






NVSS J2135 10-333703* 


2.518±0.001 


2.1x102^ 


-1.45 


47 ± 2 


22.5 ± 0.8 


10.7 ± 0.5 


5.6 ± 0.3 


4.1 ± 0.3 






NVSS J2257 19-343954 


0.726±0.001 


1.3xl026 


-1.62 


85.1 ± 4 


37.5 ± 1.2 


15.3 ± 1.0 


3.7 ±0.2 


2.5 ± 0.2 


1.6 ± 0.2 


0.7 ± 0.2 


NVSS J230035-363410* 


4.0* 


7.0xl02^* 


-1.66 


35.3 ± 2 


15.2 ± 0.7 


9.1 ± 0.5 


3.8 ± 0.3 


2.7 ± 0.2 






NVSS J230123-364656* 


3.220±0.002 


4.6xl027 


-1.57 


44.3 ± 2 


20.0 ± 0.8 


8.5 ± 0.5 


4.1 ± 0.5 


2.7 ±0.5 






NVSS J230527-360534* 


> 4* 


>9.2xl02'^t 


-1.38 


63.1 ± 2 


31.3 ± 1.0 


16.0 ± 2.0 


5.8 ± 0.2 


3.8 ± 0.2 






NVSS J230954-365653* * 


> 4* 


3.6x1027* 


-1.09 


33.5 ± 2 


19.3 ± 1.5 




5.8 ± 0.3 


3.8 ± 0.3 






NVSS J231144-362215** 


2.531 ±0.002 


1.3x102'^ 


-1.23 


33.7 ± 2 


18.1 ± 0.7 


10.9 ± 1.0 


5.5 ± 0.8 


3.5 ± 0.3 






NVSS J231317-352133** 


1.2* 


1.8xl028* 


-1.23 


30.8 ± 1 


16.5 ± 0.7 


8.8 ± 0.7 


4.0 ± 0.2 


3.2 ± 0.2 






NVSS J23 1338-362708 


1.838 ±0.002 


6.7x102'' 


-1.54 


36.3 ± 1 


16.6 ± 0.7 


7.6 ± 0.3 


2.4 ± 0.2 


2.0 ± 0.3 


1.0 ± 0.1 




NVSS J23 1402-372925* * 


3.450±0.005 


2.5xl02* 


-1.22 


241.3 ± 12 


129.9 ± 3.9 


71 ± 4 


38.0 ± 2.0 


27.0 ± 2.0 






NVSS J23 1727-352606* * 


3.874±0.002 


1.5xl02** 


-1.19 


108.5 ± 5 


59.2 ± 1.8 


31.4 ± 2.0 


12.4 ± 0.6 


8.3 ± 0.4 






NVSS J232058-365 157* 


1.3* 


7.2xl026* 


-1.28 


98.5 ±5 


51.4 ± 1.6 


34.5 ± 2.0 


16.5 ±0.8 


13.3 ±0.7 


9.3 ± 0.6 


3.5 ± 1.2 


NVSS J232100-360223* 


3.320±0.005 


4.8xl027 


-1.65 


34.8 ± 2 


15.1 ±0.7 


6.6 ±0.5 


3.2 ±0.3 


1.9 ±0.3 






NVSS J232219-355816* 


>4t 


>1.4xl028t 


-1.84 


58.6 ±3 


23.1 ±0.8 


14.0 ± 1.0 


5.8 ±0.3 


4.4 ±0.3 






NVSS J232408-353547* 


0.201 1±0.0004 


2.0xl02* 


-0.99 


27 ± 2 


16.3 ± 1.2 


8.5 ± 1.5 


6.7 ±0.8 


5.2 ±1.0 


5.1 ± 0.6 


2.2 ±0.7 


NVSS J232602-350321* 


0.293±0.001 


4.5xl024 


-1.07 


28.9 ±2 


16.8 ± 1.2 




7.3 ±0.5 


6.6 ±0.7 


5.8 ± 0.6 


2.0 ±0.4 


NVSS J23265 1-370909* 


2.357±0.003 


2.8x102^7 


-1.35 


64.3 ± 3 


32.5 ± 1.1 


15.6 ±0.8 


6.8 ±0.3 


4.2 ±0.3 






NVSS J234 137-342230* 


>4t 


>4.6xl027t 


-1.31 


33.9 ±1 


17.4 ± 1.0 


6.8 ± 1.0 


3.6 ±0.3 


2.5 ±0.3 






NVSS J234 145-350624* 


0.644±0.001 


3.7x1027 


-1.10 


3190 ± 160 


1823.2 ±54.7 


1062 ± 53 


490 ± 25 


353 ± 18 


242 ± 12 


82 ± 16 


NVSS J234904-362451* 


1.520±0.003 


9.3xl02'3 


-1.23 


86.6 ±3 


46.5 ± 1.8 


28.2 ± 1.4 


17.1 ±0.9 


13.7 ±0.7 


9.4 ± 1.0 


4.3 ±0.9 


NVSS J235 137-362632** 


>4t 


>8.1xl027t 


-1.19 


68.5 ±3 


37.4 ± 1.5 


18.8 ± 1.5 


8.7 ±0.5 


5.8 ±0.3 







fRedshifts are estimated using i4'(8") = 4.633 logjg (z) + 17.066 where ii'(8") is the K-band magnitude defined in an 8" aperture iDe Breuck et all2002h . Since the K — z relation is not known 
to hold beyond z A, the six sources where the relation predicts redshifts above this are listed as 2 > 4. 
fThese sources no longer meet the original USS sample threshold of a < —1.3 (see i|3.U . 
*Sources also observed in 6A array at 4.8 GHz and 6.2 GHz. 
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Table 3. Size constraints derived from the high resolution 4.8 GHz and 
6.2 GHz observations of the 16 highest redshift candidates (see J3j on 2004 
December 10. Projected linear sizes are given for sources with spectro- 
scopic redshifts from paper 11. Positive position angles of the observed radio 
axis are measured East from North. 



Source LAS PA type Linear Sizet 

" ° kpc 



Table 4. Functional form for the rest-frame SEDs fitted by log( J^j'^ ) = 
b -f alog( + /3(logQ^)^, and its associated goodness-of-fit. 



Source 




b 


a 


/3 




NVSS J001339- 


-322445* 


2.32 


-0.97 





6.4 


NVSS J002131- 


-342225 


1.51 


-1.59 


0.57 


4.1 


NVSS J002219- 


-360728 


1.42 


-1.33 


0.32 


12.4 


NVSS J002402- 


-325253 


1.79 


-1.01 





10.4 


NVSS J002627- 


-323653* 


1.85 


-1.23 





12.0 


NVSS J011606- 


-331241 


1.60 


-1.44 





2.4 


NVSS J012904- 


-324815 


1.84 


-1.37 





1.9 


NVSS J015232- 


-333952* 


2.49 


-0.82 





8.4 


NVSS J015324- 


-334117* 


1.44 


-0.81 





14.6 


NVSS J015544- 


-330633* 


1.70 


-0.91 





1.8 


NVSS J021308- 


-322338 


1.56 


-0.90 





3.8 


NVSS J030639- 


-330432 


1.89 


-1.63 





2.0 


NVSS J202026- 


-372823* 


1.81 


-1.16 





4.6 


NVSS J202140- 


-373942* 


1.44 


-0.92 





7.1 


NVSS J202945- 


-344812* 


1.85 


-0.94 





13.5 


NVSS J204420- 


-334948 


1.89 


-1.50 





1.5 


NVSS J213510- 


-333703 


2.00 


-2.07 


0.47 


2.3 


NVSS J225719- 


-343954 


2.00 


-1.77 





11.1 


NVSS J230035- 


-363410 


1.64 


-1.29 





13.3 


NVSS J230123- 


-364656 


1.81 


-1.47 





7.3 


NVSS J230527- 


-360534 


1.97 


-1.39 





1.3 


NVSS J230954- 


-365653* 


1.47 


-1.04 





2.8 


NVSS J231144- 


-362215* 


1.49 


-1.10 





2.1 


NVSS J231317- 


-352133* 


1.45 


-1.15 





2.1 


NVSS J231338- 


-362708 


1.69 


-1.53 





1.8 


NVSS J231402- 


-372925* 


2.32 


-1.06 





4.6 


NVSS J231727- 


-352606* 


2.15 


-1.28 





2.5 


IM Vv30 JZJZUJo- 


-jDj 1 j / 


1 87 
1 .o / 


— u.y / 


n 
u 


Q 9 


NVSS J232100- 


-360223 


1.67 


-1.42 





7.8 


NVSS J232219- 


-355816 


2.45 


-2.63 


0.65 


9.7 


NVSS J232408- 


-353547* 


1.32 


-0.75 





6.3 


NVSS J232602- 


-350321* 


1.37 


-0.76 





10.0 


NVSS J232651- 


-370909 


1.87 


-1.32 





3.9 


NVSS J234137- 


-342230 


1.64 


-1.30 





3.1 


NVSS J234145- 


-350624* 


3.45 


-1.11 





3.3 


NVSS J234904- 


-362451* 


1.79 


-0.90 





13.2 


NVSS J235137- 


-362632* 


1.91 


-1.22 





1.7 



*These sources do not meet the original USS sample threshold of 
OL < -1.3 (see jxTT . 
fThe degrees of freedom for the values can be obtained directly from 
Table 2. 



6 DISCUSSION 

A remarkable feature of the 37 SEDs in our sample is that they can 
be described, in the majority (33 out of 37) of cases, by a single 
power law. For the four (11%) spectra which do show curvature, 
these flatten rather than steepen toward higher frequencies. We 
note that these four flattening sources all have a^fg' < —1.3. This 
latter behaviour is indicative of either (i) the emergence of flatter 
spectrum core or hotspot components becoming important contrib- 
utors to the total flux density at high frequencies (currently, we do 
not have observations of these four sources at high enough angular 
resolution to test this), or (ii) preferential biasing toward sources 
with alif > -1.3 which have scattered into the SUMSS-NVSS 
USS sample (see §5.2 in paper I). Not one of the 37 objects in our 
sample shows evidence of steepening toward higher frequencies. 
We remind the reader that our sample was selected between 
843 MHz and 1.4 GHz in the observed frame. The corresponding 



NVSS J002402- 


-325253 


< 1 




unresolved 


< 9 


NVSS J204420- 


-334948 


< 2 




unresolved 




NVSS J213510- 


-333703 


< 1.5 




unresolved 


< 12 


NVSS J230035- 


-363410 


< 1.5 




unresolved 




NVSS J230123- 


-364656 


< 1 




unresolved 


< 8 


NVSS J230527- 


-360534 


< 1 




unresolved 




NVSS J230954- 


-365653 


55.2* 


25* 


resolved 




NVSS J231144- 


-362215 


15.4 


64 


resolved 


125 


NVSS J231317- 


-352133 


< 1.5 




unresolved 




NVSS J231402- 


-372925 


4 


61 


resolved 


28 


NVSS J231727- 


-352606 


4 


-42 


resolved 


26 


NVSS J232100- 


-360223 


< 1 




unresolved 


< 8 


NVSS J232219- 


-355816 


< 1 




unresolved 




NVSS J232651- 


-370909 


7 


13 


extended 


58 


NVSS J234137- 


-342230 


< 2 




unresolved 




NVSS J235137- 


-362632 


6 


10 


extended 





*From paper 1. 

^Linear sizes are given for those sources with spectroscopic redshifts. 



5.1 Host galaxy identification and radio structure 

Freak flooding occurred at the ATCA during our high spatial res- 
olution observations on 2004 December 10. As a result, the atmo- 
spheric water vapour content sky-rocketed and the phase stability 
at 4.8 GHz and 6.2 GHz became extremely poor. Consequently, our 
observations on December 1 1 and 12 have not been included in the 
analysis as we could not salvage the data. Unfortunately, the signal- 
to-noise from the December 10 images is not adequate to determine 
reliable polarisation properties or accurate flux densities. However, 
they are of sufficient quality to measure the angular size and fur- 
ther constrain the positions of the host galaxies which have so far 
eluded identification (paper II). Table|3|lists the angular and linear 
sizes (for sources with spectroscopic redshifts) derived from these 
observations, or a limit if the sources remained unresolved. An as- 
terisk (*) in column 1 of Table 2 shows the subset of sources in 
Tablelll 



5.2 Spectral Energy Distributions 

For each source in the sample, we derived a functional form for 
its rest-frame SED by using a 1/ct^ weighted least squares fit to 
log(5'i,/(l + z)) versus log(!/(l -|- z)). We began by fitting each 
of the data sets with first-order polynomials. The goodness-of-fit of 
each model was determined using the test. Models were rejected 
at the 2% level. The SEDs of 33 out of 37 sources were well char- 
acterised by a single power law. First-order polynomial fits were 
rejected for four sources in the sample. For each of these, we fitted 
the data sets with second-order polynomials and repeated the statis- 
tical tests. The results of the fits are given in TableHand the SEDs 
are shown in FigureQ 
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Figure 1. Rest-frame spectral energy distributions for the 37 SUMSS-NVSS USS radio galaxies in the sample, shown in order of increasing Right Ascension. 
The flux densities listed in Table 2 are plotted along with the functional form for the SED as given in Table R1 Source redshifts are labelled in the bottom left 
comer of each plot. 
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Figure 1. - continued 
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Figure 1. - continued 
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Figure 1. - continued 



rest-frame frequencies range from ~ 4 GHz for z > 3 to ~ 1 GHz 
at the lowest redshifts. We interpret this result as evidence that, 
regardless of its redshift, the SUMSS-NVSS USS sample cannot 
be described by SEDs which steepen at rest frequencies beyond 
about 1 GHz. Lower frequency observations of the highest 
redshift galaxies are, of course, required to determine the actual 
extent of any curvature between 1 GHz and 4 GHz in the rest frame. 



selected at 1400 MHz with q?Joo < -1.3 (this paper). The under- 
lying reason for these different outcomes is unclear: possibilities 
include differences in flux density limit, finding frequency and the 
spectral index threshold which is applied to the samples. We have 
recently begun investigating these possibilities and will report our 
findings in a future paper. 



This result is in striking contrast to complete samples such 
as 3C, 6C and 7C where approximately 70% of the SEDs steepen 
toward higher frequencies I Blundell et al. 1999). The fraction 
seems to be dramatically reduced in steep-spectrum selected 
samples: down to 29% for a sample selected at 151 MHz with 
a\llo < -0.981 (Blundell etal. 1998J , and 0% for a sample 



7 THE REDSHIFT - SPECTRAL INDEX CORRELATION 

As mentioned previously the mechanisms invoked to explain 
the apparent correlation between redshift and observed radio spec- 
tral index are (i) k-corrections of steepening radio SEDs, (ii) en- 
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hanced spectral aging due to increased inverse Compton losses 
against the CMB at high redshift, and (iii) an intrinsic correlation 
between low frequency spectral index and radio luminosity coupled 
to a Malmquist bias. We now consider each of these possibilities as 
explanations of the observed z — aftoo correlation which has led 
to the discovery of a large fraction of z > 3 SUMSS-NVSS USS 
radio galaxies (paper II). 



7.1 A k-correction 

If a k-correction is mainly responsible for the large surface density 
of SUMSS-NVSS USS z > 3 radio galaxies, then we expect those 
same sources to be characterised by SEDs which steepen toward 
higher frequencies. Clea rly, our results are in consistent with 
this requirement. Indeed, Chambe rs et al ] too*) also dismissed 
a k-correction for the steep spectral index of the well-known 
z — 3.8 radio galaxy 4C-I-4I.I7, after finding an essentially 
straight spectrum between 26 MHz and 10 GHz. In Figure |2| we 
show the spectral indices measured at 2.3 GHz in the observed 
frame (upper panel) and the rest frame (lower panel) for the 28 
radio galaxies in this sample with spectroscopically confirmed 
redshifts. As expected, there is very little difference between the 
two distributions, with the only exception being a steepening in the 
rest frame for three sources in the sample with flattening SEDs. 

Whilst many authors have shown that a k- correction is a 
credible explanation for the z — a corre lation iGooal-Krishnal 
Il988t lLacvetalj|T993l ICarilfi et alJIiggg h. some of these same 
studies also show that an apparent — yet weaker — correlation 
still exists bet ween redshift and rest-frame (k-corrected) radio 
spectr al index JUacv et alj[l993L ICarilli et alJll999l:lBlundell et all 
Il999l) . These samples were all selected at frequencies in the range 
38 MHz< u <408 MHz, at least a factor of three lower in finding 
frequency than the SUMSS-NVSS USS sample, and at least an 
order of magnitude brighter in radio luminosity. The differences 
in selection frequency and luminosity could explain our contrary 
result. However, at the same flux density limit, the surface density 
of SUMSS-NVSS USS radio galax ies at 2 > 3 is similar to the 
WENSS-NVSS USS sample of De Breuck et alJ Jioom which was 
selected between 325 MHz and 1.4 GHz. Since a k-correction 
does not apply to our sample, there must be another mechanism 
contributing to a z — a correlation. 



7.2 Inverse Compton and synchrotron losses 

Synchrotron physics predicts equivalent time-independent spectral 
changes due to IC and synchrotron losses. This means that cooling 
of the electron energy distribution due to both mechanisms 
manifests itself in the same way on the SED (see 0. Enhanced 
IC losses simply shorten the timescale over which this cooling 
occurs (Eguation llOL So, for example, a radio lobe with a given 
dynamical age will have suffered more radiative losses at high 
redshift due to the (1 + z)* dependence. The critical issue here 
is that this will simply lead to a lower break frequency. It will not 
steepen the SED beyond the canonical Aa — 0.5, for a continuous 
injection model. Compared with a low-redshift radio galaxy which 
is sampled beyond its break frequency, enhanced IC losses will not 
lead to steeper radio spectral indices in high-redshift sources. 

We also note that if equipartition between the magnetic field 



and relativistic electrons is a valid approximation at high redshift, 
then synchrotron losses will continue to dominate over IC losses 
out to z ^ 3 in radio lobes (for ty pical lobe magnetic field strengths 
of 35/iG; Hardcastle e t al. 2002) and out to z ~ 8 in hotspots (for 
typical hotspot strengths of 150/iG: .Hardcastle et al.2002') . 

7.3 Luminosity dependence 

Because of Malmquist bias, any mechanism which produces a 
correlation between spectral index and radio luminosity will give 
rise to an indirect cor relation between spectral index and redshift. 
iBlundell et alJ (l999) find that radio spectral indices (measured 
at rest-frame 151 MHz) are indeed steeper in more powerful 
radio galaxies. They posit that the correlation arises because radio 
galaxies with larger luminosities also have larger jet powers which, 
in turn, produce larger magnetic fields when the jet kinetic energies 
thermalise at the hotspots. The larger magnetic fields result in 
more rapid electron cooling times (Equations |4| and |9} which lead 
to steeper electron energy distributions injected into the lobes. The 
extent to which the spectral index - radio luminosity correlation 
contributes to the large surface density of USS HzRGs in the 
SUMSS-NVSS sample remains unclear; in order to explore this 
further, the degeneracy between luminosity and redshift in this 
sa mple must be broken, as it was for the complete samples studied 
bv lBIundelIetaIJ^I999^ . 



8 WHAT CAN WE LEARN FROM NEARBY CLUSTER 
GALAXIES? 

The spectral index distribution fo r a sample of nearby z < 0.5 
radio galaxies selected at 5 GHz <Kuehr et alJll98lt IStickel et alJ 
Il994h is shown in Figure |3| This nearby sample has been selected 
at a similar rest-frame finding frequency to the z > 3 radio 
galaxies in the SUMSS-NVSS USS sample. The average spectral 
index of the nearby sample is a| 5 — —0.85 with less than 
1% having indices as steep as —1.3 (the USS classification). If 
the rest-frame spectral index distributions of nearby (and less 
powerful) and distant radio galaxies are similar, this clearly implies 
that by using a spectral index cull as severe as —1.3, the high 
redshift radio galaxies being selected are not 'typical' of the 5 GHz 
selected samples, but fall at the extreme edge of the spectral index 
distribution for radio sources. 

The homogeneity of rest-frame SEDs for the 37 sources in the 
sample might imply that the physics leading to their steep spectral 
indices can also be explained by a uniform mechanism. However, 
as always the Malmquist bias inherent in single samples prevents 
us from breaking the redshift-luminosity degeneracy. Nonetheless, 
it is tempting to ask what self-consistent mechanism could at the 
same time explain the physics of USS radio galaxies, and demand 
that it correlates with increasing redshift. 

More than three decades ago it was realised that nearby 
radio sources in steep-spectrum-selected sam ples reside almost 
exclusively in r ich clusters of galaxies (Bald win & Scottlll973t 
ISlingd "' 1974a b). It is also known that cluster radio sources dis- 
play steeper spectral indices than field radio sources, with the 
steepest spectrum sources residing closest to the cluster centres 
iSleeetalj[l98l . This has been interpreted as the manifesta- 
tion of pressure-confined radio lobes which slow adiabatic expan- 
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Figure 2. Spectral index versus redshift, for the 28 galaxies in this sample with spectroscopically confirmed redshifts. Upper panel: measured at 2.3 GHz in 
the observed frame. Lower panel: measured at 2.3 GHz in the rest frame. Since all but three sources are characterised by a single power law, these plots are 
almost indistinguishable. The horizontal dashed line at a = —1.0 helps to guide the eye. 



sion losses iBaldwin & Scotll973t E ^ssarov & Gubano\il994 
Ijones & Prestonl200l h. A radio lobe will expand adiabatically un- 
til gas pressure equilibrium is reached between the lobe and the 
ambient gas pressure. From this it follows that in the centres of 
rich clusters, and other similar environments where the ambient 
baryon densities are large, radio lobes will be pressure-confined 
and lose energy primarily via synchrotron and IC losses. In more 
rarefied environments, like those surrounding isolated galaxies, the 
luminosity of the radio lobes would fade out much faster due to 
adiabatic expansion energy losses and decreasing surface bright- 
nesses. Note that these nearby sources are of the lower radio lumi- 
nosity Fanaroff-Riley type 1 (FR I) classification jFanaroff & RilevI 
1 1974.) . The high-redshift sources being selected using USS tech- 
niques, on the other hand, invariably have to have radio luminosi- 
ties comparable with with the more powerful FRII class to be 
above the flux selection limit. We note that previous studies of 
USS-selected radio galaxies fi nd a large (3 0%) fraction of compact 
steep-spectrum (CSS) sources jPe Breuck et al. 2000). In addition, 
USS radio galaxies with classical FRII morphologies at z 2 
are more distorted compared with their low-redshift FR II counter- 
parts iCarilli et al. 1997; Pentericci et al. 2000a). So, although we 
acknowledge the danger in drawing analogies between local steep- 
spectrum FRIs and distant FRIIs, we believe that high ambient 
densities play an important role in the physics leading to the un- 
usually steep radio spectral indices of both classes of source. 



8.1 A new explanation for the z — a correlation 

Nearby FRI radio galaxies with ultra-steep spectral indices are 
rare. Those which do exist reside overwhelmingly in regions of 
high baryonic densities. By analogy, we speculate that if there is 
evolution in the richness of the environments of powerful radio 
galaxies in the sense that these radio galaxies are increasingly 
more likely to occupy regions of higher ambient densities at 
higher redshifts, then we would expect to find, on average, more 
steep-spectrum radio galaxies as a function of redshift. This would 
be a natural, and physical, interpretatio n of the z — a correlation . 
We note here that almost a decade ago lAthreva & Kapahil il998h 
also posited that the z — a correlation is driven by the propagation 
of radio jets into regions of higher gas density as a function 
of redshift. These authors asserted that the injection of steeper 
electron energy spectra occurs naturally as a function of redshift 
in first-order Fermi acceleration processes due to the decreased 
hotspot advance velocities in a progressively denser (and hotter) 
inter-galactic medium. 

There are some good reasons we believe that radio galaxies 
are increasingly more likely to occupy regions of higher ambient 
densities at higher redshifts. Most importantly, cosmological ex- 
pansion of the universe means that the average gas density of the 
intergalactic medium was greater in the past by a factor of (1 + 2)"^. 
In addition, there is abundant observational evidence to support the 
notion that distant radio galaxies reside in high density environ- 
ments: 

(i) Firstly, the environments of many powerful radio galaxies 
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within the range 1 ;$ 2 < 2 show an excess of companion 
galaxies within a 1 Mpc radius iBomancini et aljl2004l : lBest et alJ 
[2003; Best 2000, and references therein). The direct detection 
of star-forming companion galaxies surrounding z > 2 radio 
galaxies, with narrow-band Lya or Ha emission line searches, 
also indicates that steep-spectrum radio galaxies live in regions 
of galaxy overdensity comparable with the richest galaxy clusters 
(Kur k et al. 2000; Pe ntericci et al. 2000; Venemans e t al. 2 002, 
l2004lMilev et al 120041: 1 Venemans et all2004l : IOverziej!200€l) . 

(ii) Secondly, radio polarimetric imaging of z > 2 radio galax- 
ies and radio-loud quasars has shown that many have rotation mea- 
sures which, if in trinsic to the source environment, are in e xcess 
of 1000 rad m"^ iCarilli et alJll997l : IPentericci et alJl2000bl) . The 
rotation measure (RM) is given by 

RM = 811.9 / UeB dl, (16) 
Jo 

where Ue is the electron density in cm~^, B is the magnetic 
field in /iG and L is the path length in kpc. Unless the magnetic 
fields are somehow unusually large at high redshift (a situation 
that seems unlikely), this implies that the extreme RMs found in 
HzRGs are the result of d ense environments sim ilar to nearby 
X-ray cooling-flow clusters JCarilli & Tavloj2002L and references 
therein). Extended X-ray emission has been attributed to thermal 
emission from dense ( 0.05 cm~^) gas in a t least one HzRG with a 
large rotation measure iCarilli et all2002h . 

(iii) Thirdly, observations of the fields surrounding high- 



redshift radio galaxies have revealed that many of these sources 
are enveloped by extended (up to 100 kpc) Lya haloes (e.g. 
jbeyetal. 2005; Weidineer et al. 2005; Reuland et al. 20ol 
IVillar-Martm et all2003f) . Others have been shown to be immersed 
in massi ve reservoirs of m olecular ga s and dust (e.g. .Stevens et all 
^003; Reu land et all2004 rPapadopou los et all2000ll200lh . These 
observations are all consistent with a scenario in which radio jets 
in high-redshift radio galaxies are propagating through denser 
environments than at low redshift. 

(iv) In addition, radio imaging of some high-redshift radio 
galax ies reveal unusual knotty radio em ission along the radio 
axis <Carilli et alJ Il997t IPentericci et alj EoOOhl This has been 
interpreted as evidence of "frustrated jets" propagating through a 
medium which is dense and clumpy on scales of between 50 and 
100 kpc. 

(v) And finally, the total (but unvirialised) mass contained 
within high-redsh ift protocluster s is estimated to be in the range 
2 — 9 X 10^'* M0 ( Venemans 2005), similar to the virialised masses 
of the richest present-day galaxy clusters. 



We have recently begun an observational campaign to test 
our hypothesis that high-redshift USS radio galaxies are located in 
environments with higher gas densities compared with their low- 
redshift FRII counterparts. This remains a work in progress. The 
results of our investigation will be given elsewhere. 
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9 CONCLUDING REMARKS 

In this paper we have presented muUi-frequency radio observations 
for a sample of 37 steep-spectrum radio galaxies. The sample was 
originally designed to find distant radio galaxies by exploiting a 
well-known correlation between redshift and radio spectral index. 
We determined rest-frame SEDs for the sources in our sample, 
which cover a broad range of redshifts from 0.1 < 2 < 4.0. We 
have found that 34 (89%) of the SEDs are well-characterised by a 
single power law, with the remaining 4(11%) sources showing ev- 
idence for spectral flattening toward higher frequencies. This result 
appears in stark contrast to complete samples of radio sources in 
which approximately 70% of SEDs progressively steepen toward 
high frequencies. We have begun investigating whether differences 
in the sample selection frequencies, flux density limits or spectral 
index thresholds are responsible for these differences. We also note 
that the Malmquist bias inherent in our sample currently prevents 
us from separating out the contribution from a liuninosity-spectral 
index correlation. 

Coupled to the rarity of USS sources in nearby GHz-selected 
samples, we are not persuaded that a k-correction plays a dominant 
role in our SUMSS-NVSS USS sample. Instead, we believe that 
USS HzRGs are intrinsically steeper than their nearby counter- 
parts. Since the surface density of a > 3 radio galaxies in our 
sample is up to three times higher than in lower-frequency-selected 
USS samples, and since the derived space density of 2: > 3 
radio galaxies in our sample is consistent with representing the 
majority of all powerful radio galaxies at high redshift (paper II), 
we speculate that a mechanism other than a k-correction or 
limiinosity-induced correlation is driving the — a correlation. 

There is a well-known trend between the spectral index and 
the ambient surroundings of nearby, low-limiinosity radio galaxies, 
in the sense that the steepest spectrum sources reside at the centre 
of the richest galaxy clusters. By analogy, the z — a correlation 
could be driven by evolution in the richness of the environment of 
radio galaxies if a higher fraction of radio galaxies reside within 
cluster-like ambient densities as a function of redshift. Further in- 
vestigations are now underway to test this hypothesis. 
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